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AN EMPIRICALLY DERIVED BASIS FOR CALCULATING THE AREA, RATE, AND DISTRIBUTION 

OF WATER-DROP IMPINGEMENT ON AIRFOILS 1 

By Norman R. Bergrctn 


SUMMARY 

Ah empirically derived basis for predicting the area, rate, and 
distribution of water-drop impingement on airfoils of arbitrary 
section is presented. The concepts involved represent an initial 
step Uncard the development of a calculation technique which is 
generally applicable in the design of thermal ice-prevention 
equipment for airplane wing and tail surfaces. It is shown 
that su fficiently accurate esitmates, for the purpose of heated- 
wing design, can be obtained by a few numerical computations, 
once the velocity distribution over the airfoil has been determined. 

The calculation technique presented is based on results of 
intensive water-drop trajectory computations for five airfoil 
cases which consisted of 15-percent-thick airfoils encompassing 
a moderate lift -coefficient range. The differential equations 
pertaining to the paths of the drops were solved by a differential 
analyzer. 

INTRODUCTION 

The design of thermal iee-prevention equipment for air- 
plane whig and tail surfaces has progressed to the point where 
the amount and distribution of heat flow can be calculated 
for specified flight and icing conditions (reference 1). This 
design procedure requires information as to the area, rate, 
and distribution of water-drop impingement on tbe leading 
edge of the airfoil section being analyzed. In the past, 
area and rate of water-drop impingement have been esti- 
mated by using a method involving the substitution of a 
circular cylinder for the airfoil leading edge, as suggested in 
references 1 and 2. This substitution method is adequate 
for design purposes for some combinations of cylinder diam- 
eter and drop size, but it 'can produce sizable errors for other 
combinations (references 1, 3, and 4). 

A second means of estimating the area and rate of water- 
drop impingement on airfoils is provided' by reference 3. 
This method is more accurate than the cylinder substitution 
method, but the calculation procedure is somewhat laborious 
and, as a result, its use is not too practicable in a complete 
design study where a large number of water-drop trajec- 
tories are usually required. 

To establish a procedure which would eliminate the 
laborious computations of water-drop trajectories in the 
des ign of wing thermal ice-prevention equipment, it became 
apparent that a large number of water-drop trajectories 

1 Su mmari zes material presented In NACA TN 2176 entitled “An Empirical Method Per- 
mlttlnir Rapid Determination of the Area, Rate, and Distribution of Water-Drop Impinge- 
ment on an Airfoil of Arbitrary Section at Subeonic Speeds,” by Norman R. Bergmn, 1SS1. 


would be required for study. Experience with calculating 
trajectories by the method of reference 3 had shown that the 
pattern of water-drop impingement for drop sizes usually 
encountered in flight can be related most directly to velocity 
distribution over the surface of the airfoil. Airfoil shape 
Itself appeared to have an effect on the pattern of impinge- 
ment, but to a lesser degree than velocity distribution. Five 
airfoil cases were chosen as being the minimum which could 
be expected to provide sufficient data to include the effects 
of these two factors. Water-drop trajectories were com- 
puted for these five cases. 

This report presents some of the results of the water-drop- 
trajectory computations described in detail in reference 5 
(NACA TX 2476, 1951). In addition, the method derived 
empirically in reference 5 for rapidly estimating area, rate, 
and distribution of water-drop impingement is discussed. 
The limitations of this method and the teclmique employed 
in its use are also presented herein. 


SYMBOLS 


The following nomenclature is used throughout this report: 
a airfoil mean-line designation, fraction of chord from 
leading edge over which design load is uniform 
a* instantaneous drop-acceleration ratio, dimensionless 

A area normal to flow direction outlined by several 

trajectories at free-stream conditions, square feet 
A, area of impingement outlined on an airfoil surface 
by trajectories starting at free-stream conditions 
from an initial area of A, square feet 
c chord length of airfoil, feet 
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concentration factor 



> dimensionless 


drag coefficient of drop, dimensionless 
section lift coefficient, dimensionless 
collection efficiency of airfoil based on airfoil maxi- 
mum thickness, percent 

rate of change of velocity along the stagnation 
streamline at the stagnation point | 


dimensionless 

frontal height of airfoil, fraction of chord 
slope of airfoil contour at a particular chordwise posi- 
tion, dimensionless 
length of span, feet 

liquid-water content of icing cloud, pounds of water 
per cubic foot of air 
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weight rate of water-drop impingement per unit of 
surface area, pounds per hour, square foot 
weight rate of impingement of water drops on a body, 
per unit span, pounds per hour, foot 
ratio of the vector difference between the local air 
and drop velocities . to free-stream velocity 

|, dimensionless 

radius of drop, feet 

Reynolds number for drop at relative velocity PV 

( ? 70 

Revnolds number for drop at free-stream velocity F 

e-?) 

distance along airfoil surface from leading edge, 
positive on upper surface and negative on lower 
surface, feet 

distance along water-drop trajectory, fraction of 
chord 

time, seconds 

equivalent ellipse thickness ratio for a low-drag air- 
foil fraction of chord 



maximum tluckness of airfoil, fraction of chord 
component of local velocity parallel to chord line, 
feet per second 

local velocity of air or drop, feet per second 
component of local velocity perpendicular to chord 
lino, feet per second 
free-stream air velocity, feet per second 
rectangular coordinates for a system of axes having 
the origin at the airfoil leading edge and the x 
axis, positive toward the trailing edge, lying along 
the airfoil chord, fraction or percent of chord, 
retangular coordinates for a system of axes having the 
origin at the airfoil leading edge and the x' axis, 
positive in the free-stream direction, lying parallel 
to free-stream direction, fraction or percent of 
chord 

total airfoil-ordinate intercept established by two 
impinging trajectozies starling from infinity at a 
distance Ay. apart, fraction of chord 
distance between two trajectories at infinity, fraction 
of chord * 

distance between two trajectozies at infinity meas- 
ured in x',y' coordinates, fraction of chord 
distance between two trajectozies which start at in- 
finity and impinge tangentially on the airfoil, frac- 
tion of chord 
angle of attack, degrees 
specific weight, pounds per cubic foot 
angular displacement between local velocity and x 
axis, degrees 

kinematic viscosity of air, square feet per second 
airfoil leading-edge radius, fraction of airfoil chord 


time scale 



dimensionless 


4 ’ scale modulus ( 9 - )? dimensionless 

V 7 iTj 

\P stream function, dimensionless 

SUBSCHIPTS 

a air 

av average 

er critical 

d drop 

e effective . 

I ' lower surface 

max maximum 

o initial condition 

s condition at airfoil surface 

t tangential 

u upper surface. 

DERIVATION OF THE METHOD 

The method derived in NACA TN 2470 for calculating 
area, rate, and distribution of drop impingement assumes that 
airfoil velocity distribution is the primary factor influencing 
the patlis of water drops which approach an airfoil. This as- 
sumption is an outgrowth of experience in calculating water- 
drop trajectories by the method of reference 3, and it permits 
the study of water-drop trajectory characteristics according 
to the factoid which influence airfoil pressure distribution. 

DESCRIPTION OF PROCEDURE USED TO OBTAIN WATER-DROP 
TRAJECTORIES 

Five airfoil cases were selected as being the minimum 
number which reasonably could be expected to provide 
sufficient data, for showing the effects on water-drop trajec- 
tories of altering airfoil velocity distribution. These eases 
are listed in table A. 


TABLE A.— AIRFOIL CASES CONSIDERED IN WATER-DROP- 
TRAJECTORY STUDY OF NACA TN 2476 


Case 

• 

Airfoil 

Angle of 
attack . a 
[deg) 

Cl 

Leading- 
edge 
radius, p 
(percont 
chord) 

1 


0 

0 

. 2.07 

2 

“do_ _ _ 

a 

.32 

2.07 

3 


4 

.44 

2.07 

4 


0 

.44 

2.07 

5 


4 


l.soa 





Table A shows the systematic clianges in the variables which 
affect velocity distribution. . Cases 1, 2, and 3 were intended 
to reveal the effects of altering airfoil velocity distribution by 
changing angle of attack; case 4, compared to cases 1 and 3, 
the effects of altering velocity distribution by' the addition of a 
basic load distribution obtained by' cambering the mean line; 
and cases 3 and 5, the effects of changing general airfoil shape 
for a given angle of attack and lift coefficinet. The upper- 
and lower-surface velocity' distributions over the forwurd 
region of each of the five airfoils are shown in figure 1 . Ve- 
locity distributions for several Joukowski airfoils are used 
because the required velocity' components in the field of flow 
are mofe readily calculated than for other airfoils. It is 
noted in figure 1 that the variables selected did not result in a 
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' x, percent chord 


(41 Symmetrical 15-percent-ttilci: Jcmkowski airfoil; ci-0; a-0°. 

(b) Symmetrical 15-percmt-tliick Joo]cowski airibO; ci-0J2;u-2°. 

(e) Symmetrical IS-pensat-thipt Joukowsld airfoD; ci-A44ia-t 0 . 

<d> Cambered 15- percent-thick: Jookowskl airfoil; cj-O.M; a-0°; e— 1.0 mean lme. 

(e) N’ACA 55r-015 airfofi; ep=0.44;<x=4°. 

Figure 1. — Airlbil velocity distributions for the five alribn cases comprising the differential analyzer study. 

wide variety of velocity distributions, but it is believed that tion. The range in values of and B v resulting from a 

these distributions are representative of cases in which there combination of each minim um value and a combination of 

are no marked nose-pressure peaks. each maxim um value of the three constituent variables is 

The water-drop-trajectory computations were made to about 150 to 20,000 for $ and about 35 to 1,000 for By. 

encompass a speed range of 100 of 350 miles per hour These ranges in $ and B r encompass most possible combina- 

(assuming incompressible flow), a drop-diameter range of 20 tions of values of speed, drop size, altitude, and chord length. 

to 100 microns, and a variation in altitude from sea level to The problem of obtaining area, rate, and distribution of 
20,000 feet. Airfoil chord length was varied from 3 inches water-drop impingement on an airfoil is one of deter m i ning 

to 30 feet. These variables were combined into the dimen- the solution to a set of simultaneous differential equations 

sionless parameters, f and B r , which then were used. as the yielding the trajectory or path which a water drop will follow._ 
independent variables throughout the trajectory compute- These equations, a derivation of which may be found in 

272483—64 69 
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reference 6, are essentially those which result from imposing 
conditions of dynamic equilibrium on a drop moving in an 
air stream. In dimensionless form, the equations are 


d(Utl T 7 ) p CJt(U a u d \ 

(1) 

(ir Ry 24 \V v) 

d(vtlV) p C e R (v a vi\ 

(2) 

dr Ry 24 \V v) 

(0-(v-v)'+(?-v)’ 

(3) 


Basically, equations (1) and (2) define the acceleration of 
a drop at any instant in orthogonal ( x and y) directions. 
Consequently, a double integration of these equations, start- 
ing from a selected initial point ( x 0 , y 0 ), yields z and y 
coordinate values of a drop trajectory. Equation (3) is a 
simple identity used in the solutions of equations (1) and (2). 
In performing the integrations, knowledge of the quantity 
CtRj 24 (the ratio of the actual drag coefficient to that given 
by Stokes' law of resistance) is required; also required are 
magnitudes of the air-velocity components uJV and vJV as 
a function of drop location relative to the body. (See refer- 
ence 6.) Variation of the term C d R/24 with local Reynolds 
number if was taken from reference 7, while the variation of 
the air-velocity components uj V and vJV .throughout the 
flow field was obtained analytically for the Joukowski air- 
foils. In the case of the NACA 652-015 airfoil, however, 
the velocity distribution throughout the flow field was 
obtained by an electrolytic analogy technique. 3 

In carrying out the differential analyzer computations for 
the five airfoil cases, the general procedure was to assign 
values to the terms p and R v in equations (1), (2), and (3), 
to establish initial conditions, and then to obtain the water- 
drop-trajectory traces from the analyzer. For each com- 
bination of p and R v selected, several trajectories were 
traced until the two trajectories were found, one for the 
upper surface and one for the lower surface, which were 
tangent to the airfoil surface at the point of drop impact. 
The importance of these two tangential trajectories lies in 
the fact that all drops between the tangential trajectories 
hit the airfoil and all drops outside will miss. In some cases, 
after the tangential trajectories were established, the dis- 
tance between them was divided into sis: approximately 
equal spaces, and trajectories started at the boundaiy of 
each space were traced. These intermediate trajectories 
were used to obtain an indication of the distribution of 
water-drop impingement over the airfoil surface. 

WATER-DROP TRAJECTORY DATA 

In the water-drop-trajectory study, trajectories were 
calculated for assigned values of the independent variables 
p and Ry. These trajectories provided values of trajectory 
starting ordinates and surface positions of drop impingement 
from which values of the dependent variables, area, rate, and 

1 The technique o t eJeetrolytlo analogy Is based on the fact that the stream llnca in an 
Inviscld fnenmpressible fluid and the &qul potential lines in an electrical field are frovemed by 
the same equations. By means of this analogy and suitably constructed apparatus, velocities 
at any pofnt Jn the flow field around a body can be measured directly. 


distribution of impingement, could be tabulated. A typical 
set of trajectories is shown in figure 2, and the numerical 
results obtained for the five airfoil cases are presented in 
tables I through V. 

To obtain general trends from the water-drop-trajectory 
data, consideration was given to the desirability of develop- 
ing a method for rapidly estimating values of area, rate, and 
distribution of impingement that would require only informa- 
tion which readily is obtainable for any airfoil profile. 
Airfoil contour and velocity distribution were taken as the 
information available for use in a design study. This report 
develops fairly simple and direct linking of the dependent 
variables, area, rate, and distribution of impingement, to 
airfoil contour and velocity distribution. The sequence' in 
which aii-foil contour and velocity distribution most readily 
are related to the dependant variables is as follows: (1) area, 
(2) rate, and (3) distribution of impingement. Develop- 
ment of the generalizations will be presented in this order. 

TRENDS OBSERVED IN AREA OF WATER-DROP IMPINGEMENT DATA 

In order to determine the area of water-drop impingement 
on the leading edge of an airfoil for specified meteorological 
and flight conditions, the values of s/c for the trajectories 
which impinge tangentially on the upper and lower surfaces 
must be obtained. In computational methods like those of 
references 3, 6, and 7, the procedure essentially has been to 
select values of p and R r and then to determine the trajec- 
tory. Various trajectories are computed until the tangential 
trajectory for the upper and lower surfaces is found. The 
two tangential trajectories determine the farthest positions 
of drop impingement on the airfoil surface for the selected 
values of p and R v and permit calculating area of impinge- 
ment from the equation 



In the method derived in NACA TN 2476, the reverse 
procedure is employed; that is, a point on the airfoil is 
selected (s/c) and the corresponding p, and R v values which 
are associated with the tangential trajectories at that point 
arc determined. The nature of the relationship between s/c 
and the parameters p and R v is shown in figure 3. Data for 
the .figure are those of table IV for the cambered airfoil at 
zero angle of attack and a lift coefficient of 0.44. From 
figure 3, it can be seen that any specified value of s/c in the 
figure can correspond to an infinite number of combinations 
of the variables Ry and p. Consequently, it becomes neces- 
sary to select values of two variables and to solve for the 
third. In the derivation of the procedure for estimating 
area of impingement, values of s/c and R y are assumed and 
corresponding values of p are computed. 

If, the data of figure 3 could be made available for all 
airfoils of interest, the problem of determining s/c for various 
values of j (/ and Ry would not exist because the information 
obviously would be known. Because obtaining such data 
for all airfoils is impractical, the, problem in the general 
case arises in determining values of p for given values of 
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Kic.urji 3— Typical relation between farthest position of drop impingement, scale modulus, 
and freo-stream drop Reynolds number; 16-percent- thick cambered JoukowskI airfoil; 
cj-0.44; a—O*; fl-1.0 mean line. 


s/c and R v . To determine an expression for equations 
(1), (2), and (3) are utilized to give 



where 

Equation (4) expresses generally the relation between i p 
and R v at all points in a trajectory, and, therefore, it is 
applicable at the airfoil surface for an arbitrarily selected 
value of s/c which corresponds to some particular tangential 
trajectory. It remains to establish the values of Cy?/24, 
Rjltr, and a d for the selected value of s/c. Actually, since 
CiR/24. is a known function of R, the problem reduces to 
approximating R/R v and a t at the airfoil surface. 

Evaluation of R/B v at airfoil surface. — To determine 
R/R v the method of this report is based on a graphical solu- 
tion utilizing the hodograph plane. A typical plot in the 
hodograph plane of the data from the differential analyzer 
is shown in figure 4 for the cambered Joukowski airfoil. 
To show the general relation of drop velocities to air velocities 
the hodograph of air at the airfoil surface is also shown in 
figure 4. Hodograplis for tire five airfoil cases, of which 
figure 4 is an example, revealed that tlie velocity com- 
ponents for ah drops, regardless of the combination of f and 
R r , can be represented by one faired curve. In addition, it 
became apparent that the hodograph for the drops, for both 
upper and lower airfoil surfaces, always passes through the 
point Ut/V— cos a, Vd/V—s'm a. In the simplest case of an 
airfoil at zero angle of attack, the hodograph of the drops 
always, passes through an abscissa value of unity since the 
point corresponds physically to the point of maximum air- 
foil thickness where the tangential trajectories are straight 
lines and impinge upon the airfoil with free-stream air 
velocity. The coordinates at the origin of the air and drop 
hodographs correspond, of course, to the airfoil stagnation 
point. 



* -velocity component of air, and of drop, ^ 


Flora* 4.— Typical hodographs of tangential-trajectory velocltlci and air velocities on an 
airfoil surface; 16-poroent- thick cambered JoukowskI airfoil; ci»0.14; «-0°; a-1.0 mean line. 

To show the connection between the physical and hodo- 
graph planes, figure 5 is presented. Figure 5 (a) depicts 
several water-drop trajectories in the physical plane im- 
pinging tangentially at the same point s/c on an airfoil which 
is at an angle of attack a. For constant s/c (fig. 3) there arc 
an infinite number of particular combinations of ^ and Ry 
which are affine to any particular position of tangential drop 
impingement (s/c),. In figure 5 (a), a single vector repre- 
senting the drop velocity' for ah the trajectories is drawn 
tangentially to the airfoil at the point of drop impingement. 
Only one vector is shown because the tangential trajectory 
hodograplis, such as that presented in figure 4, indicate that 
ah drops impinging tangentially at a common point may be 
considered to have the same velocity. Also shown in figure 
5 (a) is a vector representing the air velocity at the point of 
tangency for the trajectories. The angle between the 
drop- and air-velocity vectors and the x axis is designated 
by the angle d. In figure 5 (b), a typical air and drop 
hodograph is shown and the same vectors as shown in the 
physical plane are indicated. The difference in length of 
air and drop vectors at a particular s/c position is numerically 
equal to the value of R/R v given by equation (3). This 
equality provides a basis for predicting R/Rv, and forms the 
starting point for the empirical method. 

Because an examination of the drop and air hodographs 
for the five airfoil cases showed that a single value of R/Rv 
can be considered to be associated with any particular s/c 
position,- the assumption is made that other airfoils will dis- 
play tins same characteristic. In order to calculate 1 values 
of R/Ry for an arbitrary airfoil, however, both hodograplis 
of the air and of the tangential trajectories are required. 
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Flora* 6. — Relationship between physical and hectograph planes tor drop and air veloc- 
ities at air ton surthce. 


The hodograph of the air velocity at the airfoil surface is 
easily obtained from the velocity distribution over the air- 
foil, so the problem is to determine the shape of the hodo- 
graph for the tangential trajectories. From physical con- 
siderations, it is known that the tangential-trajectory hodo- 
graph always will pass through the point uJV= 0, rJV=0 
and the point uJV= cos a, VdJV = sin a. 

With two points on the trajectory hodograph ahvaj-s 
known, it was postulated that, if one more point could be 
established, preferably where the vertical-velocity com- 
ponent reaches the maximum value, ‘the general shape of 
the trajectory hodograph might be reasonably estimated. 
It was noted from the hodograplis for the five airfoil cases 
that peak values of eJV and vJV were at nearly the same 
location on the airfoil surface; that is, values of i« Mr /T r 
and Vi ma jV seem to fall on a straight line through the origin. 
A comparison was made, for the five airfoil cases, of values 
of the vertical component of relative velocity between 
drop and air attained at the position of maximum vertical 
air velocity. For this comparison, values of (ty t i _.J T r ) — 


{Vi J F) and were obtained from the five airfoil cases 

and these are plotted in figure 6. An inspection of the data 
in figure 6 shows that the four Joukowski airfoil cases provide 
a simple relation between and v tt _„J V- 

By use of figure 6, a third point on a trajectory hodograph 
can be ascertained which in turn permits the general shape 
of the hodograph to be estimated. 

The point plotted in figure 6 for the NACA 65 a -015 
airfoil upper surface does not lie on the curve established by 
the Joukowski airfoil data, and a question® arises as to 
whether this difference is real. While this question cannot 
be resolved until further data are available, qualitatively, it 
would seem that the tangential-drop velocities should tend 
to approach more nearly the surface-air velocities in the 
case of low-drag airfoils because these shapes are not so 
conducive to altering the paths or speed of water drops. 

As an aid in discussing the construction of the drop 
hodograph using only three points, figure 7 is presented. 
In figure 7 the air hodograph is first drawn, and the point 
r„ aa5 /F is established. Then, of, the three methods con- 
sidered, one procedure to obtain a drop hodograph uses the 
maximum vertical velocity of the tangential-trajectory hodo- 
graph This value is determined as being less than 

by the amount (.Ca ma JV) — (c iwu J^ , ) in accordance 
with the curve in figure 6. The value of c dm JV so determined 
is assumed to lie on a straignt line connecting the origin 
and v amtt JV. The. position of r d ma JV along the radial line 
determines the value of at that particular 

•tax 

position. Values of R/Rv for other s/c positions might be 
taken, as a first approximation, as being in the same ratio to 
the air velocity at the particular s/c position as the value of 
RIR V at Cg_„JV is to UJV at r c _,J F (curve A in fig. 7). 
Thus, an expression for R/Ry at an s/c position would be: 


R IT (R[Rv)v. 

MX 


(5) 


Values of R/Ry calculated by equation (5) usually are too 
large near point X (fig. 7) where it is known that u<JV=c os a, 
vJV=sia a, so that a drop hodograph so constructed probably 
would not pass through this point, and it should. To over- 
come this discrepancy in the drop hodograph as computed, 

RIR 

assuming a constant value of jjjy based on the peak point 

of the air hodograph, a curve without reflex is faired tan- 
gentially into this drop hodograph from the point u 4 /F=coso', 
r ifY =sin a. The combination of the proportional curve and 
the faired curve comprises the drop hodograph, which is 
labeled curve B in figure 7. For the five airfoil cases 
maximum deviations between the drop hodographs obtained 
by the foregoing method and actual drop hodographs were 
of the order of 15 percent in the value of UJV. 

Two other methods were considered for establishing drop 
hodographs. One of these methods assumed R[Ry to main- 


• Some variation In the value Mil™ can be obtained by the choice of curve 

need for the drop hodograph. In the case of the NAC A 65j-0I5 airfoil, the latitude of choice 
tor a hodograph was fairly great because of some discrepancies In the velocity -component 
data corresponding to small values of i/e. The hodograph finally chosen, and which gives 
rise to the questioned point In figure 6 , Is based only on the most reliable velocity-component 
values torn the data. 
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Fiqwbe 6.— Variation of velocity difference between drop add air with maximum y-vefodty component of air for tbe live airfoil cases Investigated. 



Curve A based on x TfrWP~ ~ 

Rr \Rv/ r m (Wkjf. 

■ai **» 

Carve B faired tangentially to carve A from point X. 

Curve O based, on (it/Kr) -constant. 

Curve D based on constant. 

XitiUHE 7. — Illustration of three possible techniques for the construction of a drop hodograph 
from a specified air hodograph. 


tain a constant value equal to the valuo prevailing at the 
point u lt IV=cos a , vJV =sin a. The other method assumed 

Tf/D 

the ratio to maintain a constant value determined by 


the value of R/R v and UJV at the point uJV= cos a, 
v d /V=s'm a. The drop hodographs given by each of those 
two methods also are shown for the example in figure 7. The 
curves are labeled C and D, respectively. These two 
methods have the advantage of not requiring the use of the 
hodograph and figure 6; however, they arc considerably 
more inaccurate (maximum deviations from the drop hodo- 
graplis for the five airfoil cases being in the order of 30 
percent), due to the neglect of factum of apparent influence on 
the drop trajectories. Either one of these latter two methods 
might he useful for particular airfoil cases winch happen to 
fall considerably beyond the scope of the data used to obtain 
figure 6. 

After the tangential-trajectory hodograph lias been estab- 
lished in relation to the hodograph for air, values of R/Ry are 
available for various cbordwise positions on the airfoil. 
These values are used in equation (4) for arbitrarily selected 
values of R v and s/e. Once values of R y arc selected, 
values of R are ascertainable. Furthermore, the term 
CtfZ?/24-is the function of R tabulated in table VI. Thus, 
to solve equation (4), the only additional term to be 
evaluated is a d . 

Evaluation of the drop-acceleration term a d . — The remain- 
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ing term to be evaluated in equation (4) is the acceleration 
of the drop at the airfoil surface a t . To determine the 
variation of this term with chord wise position, values of a d 
were calculated from the trajectory data by equation (4) for 
each [of the airfoil cases presented in tables I through V. 
The procedure used in making the calculations was to com- 
pute the value of R,'Ry by utilizing values of the orthogonal 
drop-velocity components from tables I through V for 
corresponding values of ip and R v . The term was calculable 
through knowledge of RJRy and R v . The terms R/Ry, 
CtRjVA, R r , and ip were then substituted into equation (4) 
andjsolved for a d . The results for a typical case (15-percent- 
tliick cambered Joukowski airfoil) are presented in figure 8. 



Fu'.cs* 8.—' Typical chordwise distribution ol instantaneous drop-acceleration values, far 
tangential trajectories at instant of drop Impact; 15-percent- thick, cambered JooiowsM 
airfoil; et-0.44;oc—u‘‘;a-l.0 mean. line. 


Figure 8 exemplifies that drop acceleration at the surface 
of the airfoil, like the hodograph of drop velocities for tan- 
gentially iiqpinging trajectories, can be considered a single 
relation regardless of the combinations of ip and R r . How r 
the singular nature of the acceleration values arises can be 
shown as follows: 


Equation (4) may be written 



However,. since the term (R/R r )t is taken to be constant for 
a given position on the surface, equation. (6) may be written, 
for any given chordwise position, 

<z*= (const) \ pCt (7) 

Thus, according to equation (7), if the product of ip and C t 
remains constant for various values of R r at a given chord- 
wise position, then the value of a t also will remain constant. 
Comparisons were made, for the five airfoil cases, of ip C d 
products for given s/e positions over a wide range in ip and 
R r values. These comparisons showed that, for a given 
a/c position, the product of ip and C t generally is of similar 
magnitude. A sample of such a comparison for the 15-per- 
cent-thiek cambered Joukowski airfoil at 0° angle of attack 
is shown in table B in which values of ip, for chosen values 
of R v add s/e, were taken from curves faired from the data 
tabulated in table IV. On the basis of comparisons of ipC d 
products for the five airfoil cases, the assumption that a t 
is constant for a particular chordwise position seems fairly 
well justified. 

TABLE B. — COMPARISON OF PRODUCTS OF SCALE 
MODULUS AND DROP DRAG COEFFICIENT FOR A 16- 
PERCENT-THICK CAMBEfeED JOUKOWSKI AIRFOIL 


[a— 0"; ci— 0.44; *— 1J> mean line! 


*Ci 

Upper surface 

\sr 

15 

64 

25* 

1024 


62 

55 

65 

66 


152 

185 

166 

154 


*45 

348 

"337 

340 


877 

885 

880 

830 


3500 

SI 10 

2742 

2000 


8950 

7577 

6630 

6590 

Lower surface 

-20.0 

55 

67 

55 

mm 

-15.0 

271 

234 

204 



547 

510 

492 


-5.0 

1620 

1040 

1693 


-2.5 

5800 

5400 

5850 



After inferring that the value of a d can be considered as 
being unique at any particular chordwise position, regardless 
of the values of -p and By, the problem of evaluating drop 
acceleration becomes one of determining the appropriate 
value of a t to assign to each value of s/c. 

In approximating the drop acceleration at a point where 
the drop trajectory is tangent to the airfoil surface, several 
procedures were tested, as was the ease with the term R[R V . 
Of the various procedures investigated, the one which will 
be presented herein is considered most acceptable because 
the resultant accuracy is commensurate with that produced - 
by the most accurate procedure presented for obtaining 
R/Ry. In addition, the procedure is simple in application. 
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For this procedure, the approximation is made that the 
tangential acceleration of a drop at a given point on the sur- 
face is the same as the acceleration of the air along the airfoil 
surface at the same point. 4 The equation used to express 
the drop acceleration in terms of air velocity at the airfoil 
surface is: 


T U.„d<JJJV) 
a “~V x d(s/c) 


( 8 ) 


The velocity-gradient term in equation (8) can be evaluated 
simply by plotting UJV against s/c, and obtaining the slope 
of the curve at the desired s/c positions. 

Results typical of those obtained by using equation (8) 
to approximate values of a d are shown in figure .8 for the 
cambered Joukowski airfoil. The calculated points are 
denoted by square symbols. Figuro 8 illustrates the general 
finding that equation (8) provides over most of the airfoil 
lower surface values of a d which are in good agreement with 
the data. On the airfoil upper surface, equation (8) pro- 
vides drop-acceleration values which aro in fair agreement 
with the data near the airfoil leading edge; but farther aft, 
the ability of equation (8) to predict appropriate values 
diminishes appreciably. This decrease in accuracy was 
most pronounced for the Joukowski and NACA 65 2 -015 
airfoils at 4° angle of attack. For the two 4° angle-of- 
attack cases, the inability of equation (8) to represent actual 
drop acceleration values- fairly far aft on the airfoil surface 
apparently is because the drops impinging in this region 
have, sufficiently large inertia so as not to respond to the 
very rapid changes in surface-air velocities prevailing near 
the position of maximum air velocity. Except quite near 
the leading edge, the trajectories are fairly straight, indi- 
cating that the impinging drops do not respond appreciably 
to the vertical components of air velocity. Thus, another 
approximation of drop acceleration can bo obtained by 
using the x components of air velocity. In equation (8)> 
UJV would be replaced by uJV so that 


a d =ylX 


d(uJV) 

d(s/c) 


(9) 


Results obtained by using equation (9) are presented in 
figure 8 using the cambered Joukowski airfoil as a representa- 
tive illustration. The values calculated by equation (9) aro 
shown in the figure by triangular symbols. For the airfoil 
upper surface, the agreement between calculated values and 
trajectory data i3 good fairly far aft oh the airfoil: on the 
lower surface, the agreement also appeals to be reasonably 
good. Apparently then, equation (9) can be helpful when 
estimating a d values for airfoils at angle of attack. 

The question arises as to whether it would be possible in 
the general case, when the differential analyzer data points 
shown in figure 8 were not present, to detect the inadequacy 
of equations (8) or (9) to represent the correct values of a d . 
In this regard, it should be noted that s/c values for a„=0 

4 Only tb« tangential component at drop acceleration needs to be approximated since the 
normal component ol drop acceleration Is equal to soro at the point of tangency. That the 
normal acceleration of the drop is xero at this point can be shown by writing the equations 
expressing dynamic equilibrium of a drop. The terms involving the drop and air velocities 
are resolved normally and tangentially. A substitution of the boundary conditions at this 
point shows that the normal acceleration most equal taro. 


always can be selected because these values correspond to 
ehordwise positions of tangentially impinging straight-line 
trajectories having maximum s/c intercept. Those particular 
trajectories always can be established by constructing lines 
tangent to the upper and lower surfaces of the airfoil parallel 
to the froo-stream direction. With s/c values for a d = 0 
established, there would be some indication of when these 
equations could not truly represent the correct curve. Be- 
cause, for an arbitrary airfoil case, there is no absolute 
assurance that either equation (8) or equation (9) will 
provide values of a d which will represent the correct curve, 
it is suggested that both equations be employed in esti- 
mating values. If, in using equations (8) and (9), the value 
of s/c for which a t = 0 is found to differ materially from the 
value given by straight-line trajectories impinging tangen- 
tially on the airfoil, then the calculated values should bo 
regarded with some skepticism. In such an event, reliance 
should be placed mostly on the values of a d calculated by 
equation (8) for small s/c values, and a curve faired from 
these values to a value of zero acceleration at the known 
extreme position of drop impingement. 

Calculation of scale modulus f for s/c at the stagnation 
point. — The two preceding subsections have presented 
approximate methods by means of which equation (4) can 
be evaluated to obtain values of ^ for selected R y values at 
chosen positions on the airfoil surface. However, a special 
procedure for evaluating ^ at the stagnation point is neces- 
sary, since equation (4) cannot be used to evaluate the 
scale modulus at or very near the stagnation point. This 
procedure is more suitably discussed in connection with the 
section on rate of impingement .which follows: 

TRENDS OBSERVED IN RATE-OF-IMPINGEMENT DATA 

Another quantity of interest to tlie designer of an aircraft 
tliermal-ice-prevention system is weight rate of drop impinge- 
ment on an airfoil. An expression for weight rate of drop 
impingement per unit length of span, according to reference 
8, is given by 

Af,=3600 VmAy 0 \ (10) 

In order to evaluate the rate of impingement M, in ac- 
cordance with equation (10), the term Ay 0 ' t must he known. 
When methods like those of references 3, 6, and 7 are em- 
ployed, A y 0 ' t can be determined directly from the calculated 
trajectories which impinge tangentially upon the airfoil. 
For a procedure in which trajectories themselves are not 
determined, however, evaluation of A y/ must be based 
upon quantities which are known. 

Evaluation of A y 0 ' ( using airfoil ordinates as an inter- 
mediate parameter. — Preceding sections have shown that 
(a/c) B< and (s/e) {j can be established as a function of for 
various values of J? v ; hence, the airfoil ordinates correspond- 
ing to the farthest position of drop impingement on the upper 
and lower surfaces y Ul and y (f also can be ascertained as a 
function of ■// for various values of R v - Because values of 
y„ t and y [( can be obtained readily for a wide range of ^ and 
R v values, the data were examined for a relationship involv- 
ing Ay Ct (for small angles of attack, A y„ t is approximately 
equal to y„ t ,) and the quantity y u —yi t which will he called 
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Ay t . In this regard, Ay, t was compared with Ay, for the 
values of ^ and R y values presented in tables I through V 
for the five airfoil cases. Results typical of the comparisons 
for the five airfoil cases are shown in figure 9 for the 15- 
percent-thick cambered Joukowski airfoil at 0° angle of 
attack: 



Fio.l-SK 9.— Typical variation of the ratio of trajectory starting ordinate difference to total 
airfoil ordinate Intercept as a function of scale modulus and free-stream Reynolds number; 
l fi- percent-thick cambered Joukowski airfoil; ej-0.44; 0°; a -1.0 mean line* 


.In inspection of data for the five cases showed that the 
ratio of Ay 0( to Ay, can be considered linear with respect to 
the log of the scale modulus $ for various R v values. The 
linearity was found to exist for values of (Ay „/ Ay), <0.8 for 
the Joukowski airfoils, and for values of (Ay „/ Ay), <0.9 for 
the NACA 65 2 -Q15 airfoil; but this linearity appears to he 
characteristic only of airfoils since cylinder data from .refer- 
ence 7, when plotted in the same maimer do not show tins 
property. Of special interest in figure 9, however, is the 
fact that the ratio (Ay„/Ay), must become zero at some 
particular value of for a given value of Ry. This “critical” 
value of can be calculated from an aerodynamic property 
of the airfoil. According to references 7 and 9, for sym- 
metrical bodies at 0° angle of attack, the critical value of f 
(i. e., the maximum value for a given value of R v for which 
drops just impinge on the body) is given by 

* r »4* r ^| (U) 

OX [¥.o 


For symmetrical bodies at an attitude other than 0°, or for 
unsymmetrical bodies at an arbitrary attitude, the same 
form of equation (11) applies, but with the notation slightly 
altered ; thus, 


'Pcr — 'iRy 


WalV) 

dS 


*-0 


( 12 ) 


This change is made because the small drop which impinges 
only at the stagnation point of the airfoil follows the stag- 
nation streamline which, in the general case, is not a line 
parallel to the airfoil chord line. For simplicity, equation 
(12) shall be written 

lkr=4 R r G (13) 


In order to use equation (13), the problem of assigning a 
value of G presents itself for the case of an arbitrary airfoil. 
Since the quantities s/c and E are affected only in a minor 
way by variations in (?,* it was believed that for determining 

i Calculations have shown that negligible changes In i/c and E occur for a change In ff as' 
large as 10 percent. 


G the airfoil could be replaced by a shape more amenable to 
calculation. The assumption was made that a symmetrical 
Joukowski airfoil would he representative of that type section 
having max imum thickness fairly well forward (conventional 
airfoils), and aQ ellipse representative of that type section 
having maximum thickness well aft (low-drag airfoils). 
Since the major factors influencing the value of G are thick- 
ness and angle of attack, calculations of G were made for 
symmetrical Joukowski airfoils and ellipses of different thick- 
ness-chord ratios at various lift coefficients.® The results of 
these calculations are presented in figure 10. The data in 
figure 10 (a) are intended for use with airfoils resembling 
Joukowski airfoils and may be used directly. The data in 



Fiauu 10.— V elocity gradient along the stagnation streamline at the stagnation point, as a 
function of lift coefficient and thickness ratio for two profiles. 


* No account is taken of the effect of a cambered profile on the velocity gradient G. The 
reason for neglecting this effect is that tests using an electrolytic analogy have shown that the 
effects of camber are very small In comparison with the effects of thickness, and calculations 
have shown that only large variations In G are important In affecting the values of t!c and E 


‘ , 7°4S3 — 54 70 
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figure 10 (b) are intended for use vvitli low-drag profiles; 
however, it is first necessary to establish an “equivalent 
ellipse” thickness ratio for the low-drag section being used. 
An equivalent ellipse is defined for the purposes of figure 
10 (b) as an ellipse having its leading-edge, radius equal to 
the leading-edge radius of the airfoil, and a thickness equal 
to the airfoil maximum thickness. The major axis of the 
ellipse is thus established and the ellipse thickness ratio can 
be computed. An equation expressing the thickness ratio 
of the equivalent ellipse in terms of the airfoil leading-edge 
radius and thickness ratio is: 

-(H) 

l max 

With the aid of figure 10, the value of p er for airfoils can 
be estimated for any B v value in accordance with equation 
(13). Not only does this value correspond to the condition 
of zero rate of impingement, but it also corresponds to the 
condition of zero area of impingement. Hence, the critical 
value of p can be used for obtaining an additional point for 
area-of-impingement computations, and this value will cor- 
respond to the s/c value at the stagnation point. 

While the condition of no drops impinging on the airfoil 
surface yields one point on the curves, ( AyjAy ), versus log p, 
at least one more point is required for each value of B v in 
order to establish the linear relationships as observed in 
figure 9. To locate a second point on an isopleth of B v , it,, 
is desirable to determine a value of p corresponding to a 
chosen value of (AyjAy), somewhat less than unity. The 
reason for this specification is to procure a spread in the 
values of (AyjAy), used to ‘establish the linear relationships, 
between ( [AyjAy ), and log P, for isopleths of B v . 

In developing a procedure for determining what value of p 
is associated with a specified value of (AyjAy) , on an isopleth 
of B v , the data from the five airfoil cases were examined for 
values of some parameter, related to (s/c) Ml and (s/c) ){ , which 
could be used to fix the value of p. The parameter used to 
supply the necessary values was the efficiency of drop im- 
pingement E. The relationship between E and (AyjAy), is 
given by 



Equation (15) can be derived by starting from the definition 
of E in terms of the initial drop-trajectory ordinates 




h 



-(16) 


At the small angles of attack associated with most flight 
conditions, A y 0 ' in equation (16) can be replaced by Ay 0( 
so that 


Ay 0 =Eh 


(17) 


Then, in equation (17) , if the reference dimension h is replaced 
by t naz and both sides of equation (17) are divided by Ay t , 
and tlie terms rearranged, equation (15) is obtained. 

The trajectory data for the five airfoil cases provided, for 
different values of B r , relatively constant values of E corre- 


sponding to a value 7 of (Ay,/ Ay), = 0.8. These efficiency 
values were used to obtain an average efficiency value for 
each airfoil case. Then, by using equation (15), an nverago 
value of AySma could be computed for each airfoil ease by 
using the average efficiency values and a value of {AyjAy) 
0.8. The results are presented in table C. 

TABLE C.— AVERAGE VALUES OF A yjt m „ OBTAINED 
FROM EFFICIENCY DATA FOR THE FIVE AIRFOIL CASES 
AT A VALUE OF (Ay „/ Ay) ,=0.8 


Case 

Efficiency of Impingement, E 
(percent) 

Ay. 

fa- 

num- 

ber 

Rt 

Average 
value for 
each case 

10 

39 

64 

128 

256 

812 

1024 

2048 

1 


77.5 


77.0 


75.3 


78.0 

77.0 

0.00 

2 

70.0 



74.5 


78,3 


77.5 


75.8 

.98 

3 

72.0 



73.0 


72.0 


70.5 

• 

71.8 

.90 

4 

77.0 



82.0 

__ 

86.0 



82. 0 


81.7 

1.02 

5 

55.0 

— 

30.0 


50.8 

— 

55.0 

, 

83.4 

.71 

, 


The values of A yjtnua tabulated in table C exhibit some 
variation between airfoil cases, and figure 11 is presented t-o 
show this variation when Ay ,jt max is assumed to be a function 
only of angle of attack. In figure 11, the point for the 
NACA 65 a — 015 airfoil does not lie on theeurve presented for 
the Joukowaki airfoils. If the variation of Ayji^ with 
angle of attack shown in figure 11 is used, it is possible to 



Angle of attack, a, degrees 


Fin use II.— Eatlo of as a function of angle of attack for (Ag^Am) -0.8. 

determine, for a given value of B v , an approximate value of 
p at which (Ay o /Ay),=0.8. The procedure which may be 
used for determining this value of p is shown by a hypothet- 
ical example in figure 12. From curves of (s/c),, and {8je)i t 
as a function of log p for a specified value of B v (fig. 12 (a)), 
curves of y* ( and y, ( as a function of log p are established for 
the same value of B v (fig. 12 (b)). For the relation shown 
in figure 12 (b), there is a value of Ay,/t max which is the same 
as would be chosen from tho relat ion in figure 1 1 correspond- 
ing to the airfoil angle of attack. This particular value of 
Ayjtmai corresponds to the p value at which (Ay„/Ay),= 0.8 
for the particular B v value chosen (fig. 12 (c)), and the 

1 The procedure utilized was to determine from curves of (Ay, /Ay), as a function of log yf 
(flg. 9) the value of if- at which (Ay,/Ap)i-0,8 tor different values of Rv. Then, data from 
tables I through V were used to establish curves of E as a function of log g. for tho same 
values of Rv. On the efficiency curves, the value of E corresponding to (Ajr^Ag)i-DA for a 
particular value of Rv could be determined by locating, for the same Rv value, tbs value of 
which was established from curves, similar to that in figure 9, to correspond to (Ag^AtrL-O.i 
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<a> Farthest position ot drop Impingement. 

(b) Airfoil ordinate at farthest position, of drop tmplngement. 
mi Ordinate-Intercept ratio. 


Fiona* 12.— Graphical representation of the procedure used to obtain a value of f corre- 
sponding to Oi/JAiYt-OS 

second, point on an isopleth of R f for (Ay„/Ay) t as a function 
of log yp is thereby determined. 

The previous discussion has shown how values of Ay 9i 
may be obtained for various yp and Rv values. However, in 
the design of a thermal ice-protection system, by the method 
discussed in reference 1, it is sometimes more convenient to 
determine the rate of water-drop impingement by using the 
airfoil collection efficiency E rather than by using the term 
Ay,/. In such circumstances, equation (10) becomes 

If, =3600 VmEt^ 

wherein E would be given by equation (15) . When equation 
(15) is used and the angle of attack is other than zero, the 
limit efficiency value corresponding to straight-line trajec- 
tories will be greater than unity because h usually is some- 
what greater than f a „. 

TRENDS OBSERVED IN DISTRIBUTION OF IMPINGEMENT DATA 

Of secondary importance in the design of heated wings is 
distribution of water-drop impingement over the length of 
interception along the airfoil surface. Despite its lack in 
prime importance, information concerning distribution of 


water drops over an airfoil sometimes is desired and, there- 
fore, brief mention shall be piade of observations drawn 
from the differential analyzer results. 

An examination of the trajectory data did not reveal any 
direct empirical way to obtain a functional relation between 
impingement distribution, scale modulus, and free-stream- 
drop Reynolds number. It was found, however, that a 
graphical construction can be used to approximate the dis- 
tribution of drop impingement over an airfoil surface. The 
basis for the graphical procedure was found by examining 
the variation of the concentration factor * C as a function of 
s/c for various combinations of yp and R r . Two such vari- 
ations, which are typical of the five airfoil cases investi- 
gated, are presented in figure 13 for a 15-percent-thick 
cambered Joukowski airfoil at 0° angle of attack. The 



Distance along airfoil surface from leading edge, s/c, percent 


Figuix 13. — Surface distribution of water-drop Impingement tor a 16-perceut-thIck cambered 
loukowskf airfoil; er-O.-ii; n-o°; s-I.O mean line. 


curves depicting these variations in figure 13 are shown by 
solid lines. One curve is typical for combinations of 
and R v corresponding to curved trajectories, and the other 
curve is typical for the combination of yp and R r correspond- 
ing to straight-line trajectories (yp= 0, value of R v arbitrary). 
The curve for yp= 0 is obtained by drawing a number of 
straight-line trajectories to the airfoil to obtain values of 
the concentration factor' 


n dA+ 

~dA, 


08) 


and represents the locus of maximum possible values of C. 
This curve, which will be referred to as a limit curve, always 
can be obtained for a given airfoil because straight-line 
trajectories always can be reproduced, but the curve for 
values of C less than maximum cannot be obtained because 
the shape of the curved trajectories cannot be determined. 
Because of the shape of the C distribution curves noted for 
the five airfoil cases, and of which figure 13 is an example, a 
triangular distribution is considered useful in establishing a 
first approximation to an actual distribution. For a tri- 

1 The use of the concentration factor Cln the computation of heat requirement due to drop 
impingement Is discussed in reference 1. 
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angular distribution, the maximum value of C can be calcu- 
lated from the equation 

^ Eh 

" (19) 

which is developed in NACA TN 2476. The value of C max 
given by equation (19) is considered to lie on a line connecting 
the points 6'= 1.0, s/c= 0, and £7=0, and s/c for the stag- 
nation point. The values of (s/c) Bf and (s/c), t are used to- 
deiine the extremities of the triangular distribution for a 
value of 0=0. An example triangular distribution is 
shown hi figure 13 for the 15-percent- thick cambered 
Joukowski airfoil at 0° angle of attack. The distribution is 
constructed corresponding to values of ^=512 and B y =64 
and is compared in the figure to the distribution given by 
the trajectory data for the same values of \p and B v . 

The value of C max obtained from equation (19) always will 
be low. However, if tho triangular approximation is altered 
to correspond more nearly to the shape of the limit curve 
for the C values, while keeping the enclosed area the same 
as the triangular area, more accurate concentration-factor 
values can be obtained. The altering of the triangular 
distribution is an attempt to establish the locus of concen- 
tration-factor values which would be given by data for 
calculated trajectories. 

A PROCEDURE FOR CALCULATING AREA, RATE, AND 

DISTRIBUTION OF WATER-DROP IMPINGEMENT ON AN 

ARBITRARY AIRFOIL 

Previous sections have shown how trends derived from the 
water-drop trajectory data may be applied to determine 
area, rate, and distribution of impingement for an arbitrary 
airfoil in incompressible flow. The general procedure will 
now be summarized by using, as an example, the case of 
an Nx\CA 23015 airfoil at ct=0.5. 

AREA OF IMPINGEMENT 

The procedure for calculating area of impingement consists 
primarily in determining values of (s/c) U[ and (s/c) It . The 
following steps explain how the empirical relations derived 
from the trajectory data could be used to determine these 
values, and figure 14 incorporates necessary accompanying 
graphical relationships: 

Step 1. — Construct the following curves for use during the 
computation procedure: 

(a) A large-scale plot of the airfoil (fig. 14 (a)) 

(b) A plot of s/c versus x for both upper and lower surfaces 
(fig. 14 (b)) 

(c) A plot of k for various x positions (fig. 14 (c)) 

(d) Chordwise distribution of incompressible-flow air 
velocities over the airfoil surface (fig. 14 (d)). 

Step 2. — Construct an air hodograph (tig. 14 (e)) from the 
information in figures 14 (c) and 14 (d). 

Step 3. — Construct a drop hodograph (fig. 14 (f)) using as 
aids the air hodograph of step (2), fig. 6, and equation (5). 

Step 4.— Estimate values of drop acceleration at the airfoil 
surface (fig. 14 (g)) with the aid of equations (8) and (9), and 
the known condition of zero drop acceleration at the extreme 
position of tangential drop inpingement. 

Step 5. — Compute values of the scale modulus, correspond- 


ing to selected values of s/c, by using equation (4). Values 
of B/B v , a i} and Ci.B/24 employed in equation (4) arc ob- 
tained from figures 14 (f), (g), and (h), respectively. 

Step 6. — Plot curves of s/c versus f for isopleths of ll v 
(fig. 14 (i)) using the calculated points from step (5). Values 
of f for s/e=0 are obtained for this plot by using equation 
(13) in conjunction with figure 10. 



(&) Airfoil contour. 

Fir, che n.— Graphical relationships used In evaluating farthest position of Impingement for 
an NACA 2301S airfoil; ci-0.S; a-ZXP. 



(b) Variation of tic with chord wise position. 
Tic. use 14.— Continued. 



Co) Slope of airfoil contour as a function of chordtvlse position. 
Fiucee 14.— Continued. 
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The procedure for determining total rate of impingement, 
as has been explained in reference 1, consists of summing the 
rate of water-drop impingement for each of the drop sizes 
in an assumed drop-size distribution. A summation is pos- 
sible for each size of drop by use of the equation: 

Af,=3600 EV my naz 

The values of F, m, and y max are obtainable directly from a 
knowledge of the nature of the icing conditions and the 
airfoil shape. The procedure for calculating efficiency of 
impingement consists essentially of evaluating equation (15). 
The following steps, with the aid of figure 15, are intended 
to explain how the evaluation of equation (15) is performed: 
Step 1. — Establish the following relationships for use 
during the computation procedure: s/c as a function of y[c 
for both upper and lower surfaces (fig. 15 (a)), and y t as a 
function of tp for the desired values of R r (fig. 15 (b)). Figure 
15 (b) is obtained from figure 14 (i) by employing the con- 


version relation between s/e and y/c (fig. 15 (a)). In figure 
15 (b), use is made of figure 11 to establish the value of ip 
which corresponds to the value of {Ay,/Ay) t —Q.S. 

Step 2. — Construct (Ay 0 /Ay) , as a linear function of ip on 
semilogarithmic coordinate paper for the desired values of 
i?r (fig. 15 (c)). Two points are required to establish the 
function for each value of R v . One point is obtained from 
equation (13) already discussed in step (6) under area of 
impingement; the other point is obtained through the aid of 
figure 15 (b). 

Step 3. — Calculate values of impingement efficiency using 
equation (15). Values of (AyjAy), and Ay, used are obtained 
from figures 15 (b) and 15 (c), respectively. Results Of calcu- 
lations for the NACA 23015 airfoil are shown in figure 15 (d). 

« • 

DISTRIBUTION OP IMPINGEMENT 

Distribution of impingement is considered defined, as 
explained in reference 1, when values of the concentration 
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(h) as a function of R. 


FlflURE 14.— Continued, 

£ 


o 



(1) Farthest position of Implngmenent as a function of scale modulus. 
Figure 14. — Concluded. 


factor C are determined over the region of drop impingement. 
A 'summary of Hie procedure to establish these values is as 
follows: 

Step 1. — Determine a limit distribution curve of C versus 
s/c by equation (18). To evaluate equation (18), a plot of 
y c ' versus s/c is required (fig. 16 (a)) for straight-line trajecto- 
ries. Figure 16 (a) can fee established with the aid of a 
graphical construction of straight-line trajectories impinging 
on the airfoil being considered (fig. 16 (b)). A limit distribu- 
tion is shown in figure 16 (c) for the NACA 23015 airfoil. 

Step 2 . — Construct a triangular distribution of impinge- 
ment of (7 versus s/c. To establish this distribution, three 
values of C are located on the plot. One of these values is 
given by equation (19) and is located on a lino connecting the 
points 0— 1.0, a/c= 0, and (7=0, and s/c for the stagnation 
point. The other two points are located at a value of (7=0 
at values of s/c for farthest positions of impingement. Fig- 
ure 16 (c) shows a triangular distribution for the NACA 
23015 airfoil. 



C&) Variation of «/c with airfoil ordinates. 

Fiocbz 15. — Graphical relationships used in oral dating Impingement efficiency for un NACA 
23016 airfoil; ci-0.5; ot-3.0 6 . 

Step 3. — Modify the triangular distribution established in 
step 2 to conform with the general shape of the limit distribu- 
tion found in step 1. In performing the modification, the 
area contained within the new distribution curve is made 
equal to that contained within the triangular distribution. 
This condition usually results in a larger value of C nax . A 
modified distribution curve is shown in figure 16 (c) for a 
particular combination of ^ and R r . 

EVALUATION OF THE PROCEDURE DESCRIBED IN THIS 

REPORT 

The degree to which the final values of farthest position and 
efficiency of drop impingement, as estimated herein, depend 
upon the accuracy of determination of the intermediate 
quantities (R/R v )„ a d) and (7 was investigated by determining 
the effect of arbitrarily altering these three quantities a given 
percentage. By this means, the effect on farthest position 
and efficiency of impingement can be appraised for the se- 
lected changes in the three variables; also, some measure is 
obtained of the error introduced by the approximations used 
in the calculation procedure. 

When computations were made for the 1 5-percent- thick 
symmetrical Joukowski airfoil at a= 4°, and the values of 
(R/R v ) h a*, and O were altered by ±10 percent in all possible 
combinations, it was found that in no case was changing O 


Woter-drop trajatory Efficiency of impingmont, £", percent Ratio of starting ordinate difference to total Airfoil ordinate at farthest position 

starting ordinate, yj , percent airfoil ordinate intercept, (A^/A/), percent 0 f drop Impingement, y,, percent’ 
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id.) ElBriency of impingement. 

Fir.CE* 15. — Concluded, 



Distance along airfoil surface. from leading edge, s/c, percent 

1*1 Straight-line trajectory starring ordinates as a function of »/e. 

F Is rax 15.— Graphical relationships used In evaluating distribution of Impingement for an 
XAC.t 23015 airfoil; Ci-OJ;a-3.0“. 



Fmcai 16.— Continued. 



Distance tlong airfoil surface from- leading edge, s/c, percent 

(c) Distribution of Impingement. 

Figtje* 16.— Concluded. 


significant for farthest position of impingement. The com- 
bination of positive and negative changes providing the 
largest change in ^ resulted in a change in s/c of about 2-per- 
cent chord over most of the range in values of The 
approximations contributed an additional change of only 
about X-percent chord. 

For efficiency of impingement, the effect of a change in 
the term G alone was to make a change in efficiency of 
about 0.5 percent; the combination of positive and negative 
changes in {RjR v ) , and a s providing maximum change in 
$ made a change in efficiency of about 3 percent over most 
of the range in ^ values. As compared with these changes, 
the approximations led to efficiency of impingement values 
which differed from the differential analyzer values by about 
— 15 percent. 

While the foregoing values will not necessarily be repre- 
sentative for all other airfoils, they probably indicate the 
order of magnitude of error in area and efficiency of impinge- 
ment to be expected when the error in the terms (R[R r ) t , 
a*, and G can be kept within ±10 percent. Whether this 
sort of accuracy - always can be realized by using the pro- 
cedures suggested in this report can be ascertained only 
as more water-drop-trajectory data become available. 

CONCLUDING REMARKS 

Results of water-drop-trajectory data obtained from a 
differential analyzer have indicated trends which were used 
as a basis for devising a procedure for calculating area, 
rate, and distribution of water-drop impingement on airfoil 
sections of arbitrary profile. These trends are more firmly 
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established for airfoils resembling the Joukowski airfoils 
investigated than for low-drag airfoils, since the basic data 
were obtained for four Joukowski airfoil cases and only one 
low-drag section. Further water-drop-trajectory data are 
needed, particularly for thin airfoils (order of 5 percent thick) 
at liigh speeds, and airfoils at high angle of attack (in the 
neighborhood of 12°). Whether these new data would 
make it necessary to revise the concepts presented herein, 
replace, or substantiate them remains to be seen. Until 
such data are available, however, the method derived from 
these trajectory data should permit more complete and 
accurate calculations of the area, rate, and distribution of 
water-drop impingement on an arbitrary airfoil than other 
semiempirical methods. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., May 8, 1951 
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TABLE I.— RESULTS FROM DIFFERENTIAL ANALYZER 
STUDIES OF WATER-DROP IMPINGEMENT ON A 15-PER- 
CENT-THICK SYMMETRICAL JOUKOWSKI AIRFOIL 
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.004 

8 
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.997 
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32 
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.997 

. 013 

32 
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.997 
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i 

32 
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4 

32 
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* Denotes tangential trajectories 
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TABLE II.— RESULTS FROM DIFFERENTIAL ANALYZER 
STUDIES OF WATER-DROP IMPINGEMENT ON A 15- 
PERCENT-THICK SYMMETRICAL JOUKOWSKI AIRFOIL 
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i Denotes tangential trajectories. 


TABLE III.— RESULTS FROM DIFFERENTIAL ANALYZER 
' STUDIES OF WATER-DROP IMPINGEMENT ON A 15-PER- 
CENT-THICK SYMMETRICAL JOUKOWSKI AIRFOIL 
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MEElAlLH 

do 

—.042 

.5366 

—.0782 

512 

64 

—.3556 

do 

—.062 

.6465 

—.1363 

512 

64 


_do. 1 

— =112 

.8133 

—.2045 

2048 

256 

— . 3293 

Upper 1 

.015 

.4210 

.7513 

2048 

256 

—.3324 


.002 


.4553 

2048 

256 

— .3373 

Lower.. 

—.006 

.2858 

.2142 

2048 

256 

—.3432 

do 


.3543 

.0060 

2018 

256 

—.3471 

do 

—.032 

.4137 

— . 1191 

2048 

256 

—.3501 

do 

—.045 

.4956 

—.1923 

2048 

256 


_ do. 1 

—.072 

■ 6576 

—.2912 

8192 

1024 

—a 3382 

Upper 1 

.004 


.7521 

8192 

1024 


Lower 

—.005 

.1913 

.4371 

8192 

1024 

—.3441 

do. 

—.012 

.1787 

.1520 

8192 

1024 

— .3458 

dn 

—.020 

.2787 

.0669 

8152 

1024 

—.3474 

do 

—.031 

.2867 

—.2011 

8192 

1024 

— .S4S0 

. —do.' 

— v 04 i 

. 41S6 

—.3832 

1024 

16 

—.3195 

Upper 1 

.004 

.1147 

.6545 

1024 

16 

—.3544 

Lower *- 

—.035 

.3406 

—.2806 

4096 

4096 

64 

—,.3495 


.002 

.1047 

.0945 

64 

— .3530 


—.026 

.2116 

—.2157 

16384 

256 

—.3504 



.1597 

.8494 

16384 

256 

—.3536 

Lower 1 . . 

—.025 

.9746 

—.2256 


1 Denotes tangential trajectories. 
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TABLE IV— RESULTS FROM DIFFERENTIAL ANALYZER 
STUDIES OF WATER-DROP IMPINGEMENT ON A 15-PER- 
C’ENT-THICK CAMBERED JOUKOWSKI AIRFOIL 

[a— 1.0 MEAN LINE; a -0.44; cc-0°] 


ua 

J?r 

1* 

Surface 


UitY 

SilV 

4 

256 

0.0936 

Upper 1 


L003 


4 

256 

-.0586 

Lower ' 

.-*.310 

.996 



■I'k'B 

.0916 

Upper 1 

.310 


.013 


Wfcm 

-.0866 

Lower 1 

—.213 

. 999 



16 

.0866 

Upper 1 .... 



.022 


16 

%BVl . llSI 

Lower 1 

-.215 

.994 

iiisvrn^n 

8 

64 

.0818 

Upper ' 

.294 

L012 

.031 

8 

64 

.0636 

do — 

- .098 

.992 

.022 

$ 

64 

,0263 

do 


.988 


8 

64 

iBVi 

do 

. 002 

. 085 

.011 

8 

64 

-.0318 

Lower - 

-.033 

. 984 


8 

64 

-.0810 

dn I 

-.212 

.099 


32 

256 

.0775 

Upper 1 

.275 

■TfVvn 


■ 32 

286 


dp 

.092 

.989 

.042 

32 

256 


do— 

.034 

.976 


32 

256 


do 


.972 

.017 

32 

256 

-.0325 

Lower. 


.973 


32 

256 


—do. 1 

-.199 

.990 

-.016 

128 

1024 


Upper i 

.243 



128 

1024 

.0420 

do 


.976 

.075 ; 

128 

1024 


dn_ 


.954 

.054 

128 

1024 


v_ 


.943 


128 

1024 

-.0035 

Lower _ . 


.948 


128 

1024 

-.0585 

do. 1 

-.183 

.984 


1 6 

16 

.0377 


,211 

L 028 

.128 

10 

10 


Lower 1 -. - 

-.180 

.978 

yBBTTv^w 

04 

84 

.0312 

Upper . 

. -192 


.168 

fil 

64 


do._ — 

.061 

.938 

.131 

64 

04 


do 

1HI iVv 

.898 


64 

64 

-.0315 



.884 

.048 

64 

64 

-.0525 


— 031 

.891 

itfSES 

64 

64 

-.0731 

do. 1 

->157 

.971 


256 

256 

.0180 

Upper 1 

: .158 

L038 

.238 

288 

256 

.0010 

dn 


.893 

.188 

2fifl 

256 

0165 

do 


.839 

.127 

256 

266 

-.0340 

Lower — 

^BrrlrH 


.054 

256 

256 

-.0510 

do„ 


.838 

-.017 

286 

256 

-.0680 

do. . 

-.120 

.940 

-.129 

KtTV 

KTIyTh 

0 

Upper 1 ..'. — 


BWiTiH 

.357 

In 

muAm 


Lower 1 

-.086 

.883 

-.220 

128 

16 

-.028 

Upper 1 


.638 

.405 

128 

16 

-.0382 

—do .. 

u 

.727 

.288 


10 

— . Q4S0 

do 

.013 

.638 

.189 

128 

la 

-.0675 

Lower 


.649 


128 

16 


do 

-.019 

.844 

-.048 

128 

18 

-.0772 

do. 1 _ 


.838 

-.264 

512 

64 

-.0347 

Upper 1 


.611 

.487 

812 

■ 64 

-.0426 

do 

-- .026 

.834 

.350 

612 

64 

■UPJWri^M 

dO 

■liw 

.535 

1SK4 m 

512 

64 

-.0582 

Lower 


.486 

.031 

512 

64 

-.066 



.-544 

-.124 

512 

64 

-.0726 

do.' 


.754 

-.316 

2048 

256 


Upper i 

.048 

.686 

. 576 

2048 

256 

-.050 

do 

• .018 

.448 


2048 

256 

-.0548 

do 


.381 

.192 

2048 

250 

-.0594 

Lower - 


.388 

.048 

2048 

268 

-.0638 

^ dn 




2048 

256 

-.0686 

Lower 1 — 

-.038 


-.384 

8102 

1024 

-.0548 

Upper ' 

.025 

.471 

.650 

8192 

1024 

' -.0569 

_ dp. . . - 


.258 


8192 

1024 


Lower 

:jHrnM 

.178 

.082 

8192 

1024 


— -do 

1 H 

.184 

-.079 

8192 

1024 

jfcfiyo^B 

. _do 

-.012 

.211 

-.168 

8102 

1021 

-.0650 

.—-do. 1 


.484 

—.428 


16 


Upper 1 

^BBTiTT^B 

.846 

.561 


km 


Lower, i .... 


.317 

-.368 


■ 1 


Upper 

■Rif™ 

.286 

.527 


81 

-.0676 

Lower 1 — 

—.012 

. 143 

-.263 

■TStT^Bfi 

258 

-.0636 

Upper 1 


.113 

.391 

| 16384 

256 

-.0070 

Lower ' - - 


.048 , 

-.118 , 


1 Denotes tangential trajectories. 


TABLE V— RESULTS FROM DIFFERENTIAL ANALYZER 
STUDIES OF WATER-DROP IMPINGEMENT ON.AN NACA 
65j-015 AIRFOIL 

[ci-0.44; a-4°] 


* , 

Rk 

If 

Surface 

tic . . 

UijV 

eJY 

4 

256 

—a i 28 i 


0.281 

1.0023 

0.0814 

4 

256 

-.2817 

Lower 1 

-.523 

.9973 

.0704 

16 

1024 

-.1298 

Upper 1 

.267 

1.0073 


16 

1024 

-.2818 

Lower 1 - 

-.514 

.9973 

.0701 

2 

ie 

-.1396 

Upper i 

.250 

aSSum 


is 

16 

-.1646 

do 

.031 

.0881 

.0675 

Mfl 

16 


Lower 


.9847 



16 

-.2246 


-.074 

.0807 


M 

16 

-.2467 

do 


.9867 

.0753 

2 

16 

-.'2687 

dn 

-.266 

.0847 

.0742 

2 

16 

-.2909 

do.' . .. 

514 

.9927 

.0711 

8 

64 

-.1424 

Upper 1 

.240 

1.0107 

.1065 

8 

64 

1719 

do 


.0847 

.0935 

8 

64 

-.2163 

Lower 

-.050 

.9757 

.0824 

s 

64 

-,240# 

do 

-.125 

. 9807 

.0743 

8 

64 

-.2655 

do 

-.230 

.0797 

.0712 

8 

64 

-.2899 

do.' 

-.612 

.9907 

.0641 

32 

256 

1493 

Upper i 

. 309 

■WiTUFM 


32 

250 

-.1702 

do 

.023 

.9701 

.1156 

32 

256 

-.2193 

Lower 

-.052 

■EH 

.0874 

32 

256 

-.2437 

do 

-.131 

.0647 

.0743 

32 

256 

—.3685 

do 

-.249 

.9717 

^■(nvwl 

32 

256 

-.2804 

do. * — 



.0631 

128 


-. 1003 

Upper ' 




128 

BH 

-.1820 

do 

.008 

.9517 

.1375 

128 

1024 

-,2282 

Lower 

—.0C8 

.9427 

.0853 

128 

mwm 


do 

-.145 


.0042 

128 


-.2726 

do 

-.267 

.0557 


128 

BH 

-.2878 

do.' 

-.486 

.9837 

M®r7F7B 

16 

16 

1951 

Upper 1 

. 128 

1.0015 


16 

16 


Lower ■ 

-.481 

.9776 

.0383 

64 

64 

-.2030 

Upper ' 

. 092 

.9956 

.2597 

64 • 

64 

-.2136 

do 


.9166 


64 

04 

-.2635 

Lowor - 

-.048 

.8796 

.1092 

64 

64 

-.2787 

do — 

-.127 

.8075 

.0609 

64 

64 

-.3013 

do 

-.247 

.9266 


64 

64 

-.3111 

do. > 

-.417 

.0756 

.0074 

266 

266 

-.2143 

Upper' 

.047 

.0416 

.3410 

256 . 

256 

-.2343 

Lower 

prPTvr^B 

.8256 

.1964 

256 

256 

—.2543 

do 

-.040 

.8050 


256 

256 

-.2717 

do 

IllfESiS^B 

.8316 


256 

256 

-.2802 


— . 184 



266 

266 

-.2983 


-.356 

.9086 

-.0443 


1024 

-.2267 

Upper i 

.026 

.8675 

.4614 



—.2313 

rio 

.010 

.7728 

.3936 

I^BfitjtB 

K | L .7 B 

-.2492 

Lower... 

-.015 

.6716 

.1972 


B zrm 

-.2692 

do 


.6296 

.0300 

^mVTB 

Hi JT7H 

-.2782 

do 

IBBSSlB 

• .7876 

igfiaw 



-.2853 

do 

-.160 

.7496 

— 0092 


HtSH 

-.2883 

do. i - 

-.245 

.9175 

.0042 

128 

16 

—.2633 



.7020 

.1950 

128 

16 

-.2708 

Lower 

-.012 

. 0436 

.2336 

128 

10 

—.3945 

do 

-.042 

. 6395 


128 

10 

-.3028 

do. 

-.060 

.7104 


128 

16 

-.3050 

— do. 

-.104 

.7034 

-.0800 

128 

16 


do. '.. 

-.185 


— 1302 

512 

64 

-.2676 

Upper 1..., . 

.018 

.7636 

.5770 

* 512 

64 

—.2719 

do 


■ MW 

.4878 

512 

64 

-.2867 

Lower 


.4936 

.1605 

512 

■ 64 

-.2920 

__ do 

- BlUS 

.6504 


512 

64 


do 

-.070 


— 0989 

.612 

64 

ISiokLiB 

da 1 

145 

.8284 

WWfl 

2048 

256 

-.9737 

Upper' — 

.014 

.0006 

.6718 

2CHS 

256 

-.2748 

do 


.5986 

.6808 

2048 

256 

-.2781 


a 

.6216 

.1106 

2048 

256 

-.2861 

Lower 

-'• : vicrw * 

.1086 

. 2245 

304S 

-256 

-.2906 

do 

—.022 

.3906 

ISBEilB!? 

2048 

256 

-.2965 

do 

-.048 

.6486 

-.0907 

2048 

250 

—.2989 

do. 1, . 


.7584 

— . 17C8 

8192 

■ fi7| 

-.2798 

Upper 1 

.010 

.6970 

.7847 

8192 


MDFUiTH 

do 




8192 


—.2850 





8192 

BrivjB 

-.2886 


-.010 

HKmB 

.2464 

8192 

HTwtH 

-.2932 

do 


.2124 


8192 


-.2045 

do - 

IBa 

.3423 

— 0S68 

8192 

htt^tH 

— .2958 

do 

-.033 


-.1068 

8192 

HTQH 

-.2971 

do. ' 

—065 

.6765 

-.1807 


16 

-.2933 

Upper ' 


,4372 

Bfeiiil i! Bl 

mmtuTW 

16 


Lower ' 

-.040 

. 4891 

-.1007 


64 

-.2947 



.6712 

.8627 


04 



m^B 


-.1126 

16384 

256 

-.2952 

Upper 1 

.008 

.6730 

mmLaiiBI 

16384 

256 


hrnmSmm. 

— 019 

.1090 

-.0625 j 


' Denotes tangential trajectories. 
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TABLE YI— VALUES OF C d R !2i AS A FUNCTION OF R 


R 

CtRfl i 

s 

GiH/24 

0 

1.000 

500 

8.52 

.05 

1.000 

250 

7.38 

.1 

LOIS 

300 

8.26 

.2 

L 037 

350 

9.00 

.4 

L073 

400 

9.82 

.6 

1.106 

■ 500 

IL 46 

.s 

L 142 

600 

12.97 

1.0 

1. 176 

800 

15. SI 

1.2 

L 201 

1,000 

18.62 

1.4- 

1.225 

1,200 

21.3 

1.0 

L24S 

1,400 

24. 0 

1.8 

1.267 

1,600 

25.9 

10 

1.285 

1,800 

29.8 

2.5 

1.332 

2,000 

32.7 

3.0 

1.374 

2,500 

40.4 

3.5 

L 412 

3,000 

47.8 

4.0 

1.447 

3,500 

55.5 

5.0 

L 513 

4,000 

63.7 

6.0 

1.572 

5,000 

aao 

S.0 

1.678 

6.000 

95.8 

10.0 

L782 

8,000 

130.6 

12 

L901 

10,000 

166.3 

u 

2.006 

12,000 

204 

18 

2.109 

14,000 

243 

18 

2.198 

16,000 

285 

20 

2-291 

18,000 

'325 

25 

2.489 

20,000 

365 

30 

2.673 

25,000 

470 

35 

2.851 

30,000 

574 

40 

3.013 

35,000 

574 

50 

3.327 

40,000 

77* 

ao 

3.60 

50,000 

980 

so 

4.11 

60,000 

1175 

10C 

4.59 

80,000 

1552 

120 

5.01 

100,000 

1905 

140 

5.40 

L2X10* 

2234 

160 

5.76 

1. 4X105 

2549 

180 

8.15 

1-6X1 0 5 

2851 




